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Abstract 

We present a preliminary measurement of the CP asymmetry in singly Cabibbo-suppressed D'^ 
K^K^tt'^ decays and in the resonant decays —^ (j)7T~^ and K*^K^. We use a data 

sample of 79.9 fb~^ recorded by the BABAR detector. The Cabibbo-favored Df — > K~K~^it^ 
branching fraction is used as normalization in the measurements to reduce systematic uncer- 
tainties. Preliminary results of the CP asymmetries obtained are Acp{K~^ K~ tt^) = (1.4 it 
1.0(stat.)±l.l(syst.))xl0-2, Acp((/)vr=^) = (0.2±1.5(stat.)±0.8(syst.)) x 10-^, and Ac p{K^K*^) = 
(0.9 =b 1.7(stat.) lb 0.8(syst.)) x 10"'^. A preliminary determination of the branching ratio is 
lfD+ZK~Sn+) = (10-7 ± 0.1(stat.) ± 0.2(syst.)) x lO^^. 
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(a) (b) 
Figure 1: Examples of — > K~K~^'k~^ decays: (a) a tree diagram, and (b) a penguin process. 



1 INTRODUCTION 

Singly Cabibbo-suppressed (SCS) charged Z)-meson decays are predicted in the standard model 
(SM) to exhibit CP-violating charge asymmetries of the order of 10~^ T. Direct CP violation in 
SCS decays is possible in the interference between tree- level, Fig.^a), and penguin, Fig.^b), decay 
processes. Doubly Cabibbo-suppressed (DCS) and Cabibbo-favored (CF) decays are expected to 
be CP invariant in the SM because they are dominated by a single weak amplitude. Measurements 
of CP asymmetries in SCS processes greater than 0(10^^) would be strong evidence of physics 
beyond the standard model Hj. 

We define the CP asymmetry by 

o — 2 

\A\^ - A 

^cp = — (1) 



\A\' + 



A 



where A is the total decay amplitude for decays and A is the amplitude for the charge-conjugate 
decays. Acp is different from zero only if there are at least two different decay amplitudes where 
there has to be a relative weak phase and an induced phase shift between the amplitudes due 
to final-state interaction. Eq. can be expressed as a function of the branching fractions with 
the CPT invariance, T(D~^) = T(D~). Assuming that CF decays are invariant under the CP 
transformation, we use them as normalization factors in the asymmetry, 



B(D+- 


-*K+K-n+) 


B{D'- 


-*K+K--K- 


) 


B(D+- 


->K+K-iT+) 


B(D-- 


-^K+K--K- 


) 


B{D+- 


■*K+K-n+) 


B{D-- 


-^K+K-TT- 


) 



K+K--K-) 

This procedure reduces systematic errors since most of the particle identification (PID) and tracking 
errors cancel out. We also measure the CP asymmetry in the resonant decays (j)-K~^ and 

Finally, we present a preliminary measurement of the branching ratio ^^+~^^J^+^+^ . 
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2 THE BABAR DETECTOR AND DATASET 



This analysis is performed with a data sample recorded on and below the T(4S) resonance with 
the BABAR detector at the PEP-II asymmetric-energy e~^e~ storage ring. The BABAR detector is 
described in detail in Ref. |2|- The silicon vertex tracker (SVT) and the 40-layer cylindrical drift 
chamber (DCH) within a 1.5-T solenoid measure the momenta and energy deposition of charged 
particles. A ring-imaging Cherenkov detector (DIRC) is used for charged-particle identification. 
Photons are detected and electrons identified with a CsI(Tl) electromagnetic calorimeter (EMC). 
Muons are identified in the instrumented flux return (IFR), composed of resistive-plate chambers 
and layers of iron, which return the magnetic flux of the solenoid. 

We split our 89.7 fb~^ sample into a 9.8 fb~^ randomly selected subsample used to optimize 
the selection criteria and the remainder, 79.9 fb~^ sample used for the analysis. The subsample 
is used to finalize the analysis procedure including the study of systematic errors. This procedure 
limits selection bias. Furthermore, the same selection criteria are applied to the CF and SCS modes 
whenever possible to reduce systematic errors. We use 145 fb~^ equivalent of (generic) cc Monte 
Carlo (MC) events to determine efficiencies. 

3 ANALYSIS METHOD 

We reconstruct and |Sj mesons by selecting events containing at least three charged tracks. 
Tracks are required to have at least 12 measured DCH coordinates, a minimum transverse mo- 
mentum of 0.1 GeV/c, and to originate within 1.5 cm in xy (transverse to the beam) and ±10 cm 
along the z-axis (along the e~ beam) of the nominal interaction point. Kaons are identified with 
a likelihood ratio constructed with dE/dx likelihood functions from the SVT and DCH, and a 
DIRC likelihood function constructed with the Cherenkov angle and number of photons. Pions are 
identified as tracks that fail a loose kaon identification criteria. Three charged tracks are fitted 
constraining their paths to a common vertex, and accepted if the fit probability -P(x^) > !%• We 
reject and Df mesons from B decays by requiring that the D momentum in the center-of-mass 
(CM) frame be above 2.4 GeV/c. 

In order to reduce the remaining combinatorial background we construct a likelihood ratio 
(r) from the probability density functions (PDFs) of the following discriminating variables for 
the D'^ and Df mesons: CM momentum (pcm) and vertex probability (x^-based) with beam spot 
constraint (Pbs(x'^))- The signal PDFs are obtained with a background-subtraction technique from 
the data subsample. For decays, the signal band is defined as m£)+ G [1.854,1.882] GeV/cP 
and the sideband mass regions as m£,+ G {[1.819,1.833] U [1.903,1.917]} GeV/c^ [see Fig. El (g)]. 
A joint likelihood function is constructed for the signal, Csig = ni^sig(^«)! and the background, 
^bkg = ni'Cbkg(^i)' where i runs over the variables used. The ratio of the joint likelihoods r = 
>Csig/>Cbkg is a powerful variable to separate signal and background. About 16% of the events have 
more than one meson. For such events the likelihood ratio is calculated for each candidate and 
the candidate with the highest likelihood ratio is selected. 

The sensitivity S/AS, where S is the signal and AS is its error, is optimized as a function 
of the likelihood ratio. Using the subsample, the optimal selection is found to be r > 4.3. This 
criterion is applied to both CF and SCS decays. 

The resonant final states (pir^ and — > K*^K~^ are selected by requiring that the 

invariant mass of the resonant decays be within 0.01 GeV/c^ and 0.05 GeV/c^ of the nominal (p 
and K*^ masses, respectively. In addition, the signal is optimized by a selection on the cosine of 
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the heUcity angle {cos Oh). In the — > 07r+ decay mode, the helicity angle is defined as the 
angle between the K~ and the vr"*" in the cp rest frame. In the — > K*^K~^ decay mode, the 
helicity angle is defined as the angle between the K~ from the K and the in the K rest 
frame. Maximum sensitivity is obtained when |cos^i/| > 0.2 and |cos^//| > 0.3 for 
and K*^K~^, respectively. 

The CF — > K~K^tt'^ decays are selected similarly. The signal and sideband mass regions 
are chosen to be mi^x^ E [1.944,1.992] GeV/ c^, and rriKKn G {[1.914, 1.938]U[1. 938, 1.998]} GeV/c^, 
respectively [see Fig. |21 (a)]. The only difference from the SCS case is that contamination from 

— > i^~7r+7r^ decays is removed. In the KKtt candidates, the kaon with the same charge as the 
pion is labeled as a pion and then the Kirn invariant mass is calculated. We observe a peak 
indicating that part of the Df signal is composed of misidentified candidates. Events in the 
region 1.855 < mK-K-w ^ 1.883 GeV/c^ are removed from the sample. 

Contamination from D*^ — > > K~ it~^ K~ K~^)'n~^ decays is removed with a kinematic re- 
quirement on the invariant mass, rnx-h+ ^ 1-84 GeV/c^. In the case of the -D"*" K~Tr~^7r~^ de- 
cays, both Ktt combinations must satisfy the requirement. Partially reconstructed K^tt^tt^ 
events are also misidentified as events when the vr*^ is missed and the charged pion is misiden- 
tified as a kaon. Most of these events are removed by labeling a kaon track as a pion and applying 
a restriction on the mass difference 0.139 < (m^-^+^+ — mj^-^+) < 0.150 GeV/c^. 

The optimized selection criteria are applied to the final sample to obtain the signal yields. 
Figure 121 shows the mass distributions. The yields, listed in Tabled are computed by subtracting 
a scaled background estimate obtained from the sideband mass regions from the number of events 
in the signal region. This technique minimizes sensitivity to background shape assumptions. 

Table 1: Summary of yields in signal and normalization modes. 



Parent Charge + 



D+ - 


K- 


-K+7r+ 


21632 ± 228 


20940 ± 226 


D+ - 




f 


5452 ± 87 


5327 ± 86 


D+ - 


K*^K+ 


5247 ± 96 


5113 ±96 


Dt- 


K- 




23066 ± 217 


22928 ± 214 


D+ - 


K- 




236254 ± 570 


237616 ± 571 



The efficiencies needed for the Acp calculation are obtained from a sample of MC (generic) 
cc events. The selection criteria optimized for the subsample are applied to the MC sample. 
Efficiencies are then calculated as ratios of the numbers of selected signal MC events to numbers 
of generated events. The decay efficiencies are shown in Table [21 

We use Eq. (recognizing that branching fractions are proportional to yields divided by 
efficiencies) to obtain Aqp- As cross-checks, we calculate the CP asymmetries with two other 
methods (which would have larger systematic errors than the primary method): (i) using the CF 
K^TT^TT^ branching fraction as normalization, and (ii) without any normalization 

factor, Afjp, and the results are shown in Table |31 A study of the CP asymmetry in bins of the 
K~K~^TT~^ Dalitz plot indicates that the asymmetry is consistent with being constant and 

zero. 

The relative branching ratio ^Y(p +^ K -t+ ) been calculated as follows. The CF and SCS 
Dalitz plots are first binned to have equally populated bins. Then, the signal and normalization 
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yields and efficiencies are calculated bin by bin. The efficiency-corrected yields are then summed 
and divided to obtain the ratio. We obtain a branching ratio of (10.7 ± 0.1(stat.)) x 10~^ with the 
final sample. 

Table 2: Summary of the efficiencies for positively (e"*") and negatively charged I?(5)-meson 
decays. Efficiencies are in percent [%]. 





Decay 


e+ 


e 


D+ - 


K-K+n+ 


8.20±0.04 


8.26±0.04 


D+ 




7.67±0.07 


7.63±0.07 


D+ 


K*^K+ 


5.88±0.07 


5.90±0.07 


D+ 




9.90±0.02 


10.17±0.02 




K-K+TT+ 


3.77±0.02 


3.79±0.02 



Table 3: Summary of CP asymmetries measured in three different ways. 





Acp [10-2] 


[10-2] 


4% [10-2] 


{K-K+T^^) 


+1.36 ± 1.01 


+0.58 + 0.86 


+2.07 + 0.84 


(07r±) 


+0.24 ± 1.45 


-0.54 + 1.35 


+0.94 + 1.33 




+0.88 + 1.67 


+0.10 + 1.58 


+1.58 + 1.57 



4 SYSTEMATIC ERRORS AND CROSS-CHECKS 

We estimate the systematic error on the CP asymmetries in three different ways. The first approach 
combines estimates of the contributions from various identified sources listed in Table|3 The second 
and third estimates come from partly redundant direct studies of asymmetries in the normalization 
and control samples. 

The first row of Table 0] gives the error (0.06%) assigned to Acp due to small differences in 
momentum spectra of vr, K from and Df decays. This error is estimated as three times the 
maximum difference in vr, K MC-efficiency asymmetries in and D'^ decays. We evaluate an 
error for the background subtraction by changing the widths of the sideband mass regions. The 
error is taken to be the difference in the central values of Acp- The errors in the likelihood-ratio 
technique are estimated with two variants: (i) tightening the likelihood ratio to produce a 10% 
change in the decay yields, and (ii) using another likelihood ratio (ri) which incorporates a third 
variable, the distance in the j;y-plane from the interaction point to the vertex {dxy). Case (i) 
is obtained at r > 6.0 and the maximum sensitivity in case (ii) is at ri > 8.8. Table 0] summarizes 
these systematic errors for the CP asymmetries. The total errors are 0.8%, 0.7%, and 0.7% on the 

— > K~ K'^'K'^ ^ — > (/>7r"'", and K*^K~^ asymmetries, respectively. 

Our second estimate is the larger of the differences between Acp and the other two measure- 
ments A^^p and A^^^p. The error is 0.8% on both the inclusive K^K^tt^ and for the resonant 
asymmetries. 
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Our third and final estimate, which is apphcable only to the inclusive three-body final states, 
is based on CP asymmetries for the CF decays D'^ K~'k~^-k~^ and — > K~K~^'k~^ (since 
these are expected to be zero within the SM). The high-statistics K^tt^tt^ control mode is 

used to search for the scale of systematic effects intrinsic to the detector. The D'^ — > K^tt'^ 
mode, which is also our primary normalization mode, is largely insensitive to these effects since 
it has the same final state as our signal decay. In — > K~K~^tt~^ decays, both the D'^ and 
the D~ decay to two oppositely charged kaons; only pion charge differs in particle and antiparticle 
decays. In — > K~it~^-k~^ decays however, all three particles have opposite charges in particle and 
antiparticle decays. For these control samples, we obtain asymmetries of (-1-1.1 =b 0.2) x 10~^ and 
(-F0.6±0.8) X 10"2 (statistical errors only) for L»+ K' 'tt'^'k'^ and Df K K'^-k^ , respectively. 

We find that these observed control-sample asymmetries are almost entirely in the efficiencies 
derived from our simulation of the detector, rather than the signal yields (see Tables ^ and IS]). With 
corrections based on estimates of low-energy nuclear interaction effects which are not accounted for 
in the simulation we find the smaller asymmetries -1-0.8 x 10~^ and -1-0.4 x 10~^ for — > K~'k~^-k~^ 
and K'^tt^ , respectively. These results we interpret to mean that our present simulation 

of particle interactions in the material of the detector is incomplete. 

Even though our definition of Aqp [see Eq. 0] invokes a normalization chosen to eliminate 
effects that may be important here, we choose here the conservative estimate that the systematic 
error is measured by the magnitude of the departure of these CF asymmetries from the expected 
null values. 

We chose as our systematic error in the CP asymmetry the largest of all applicable estimates, 
1.1% on — > K~K~^7r~^ and 0.8% on the resonant decays. 

The CP asymmetry has been cross-checked as a function of the lab momentum as well as 
by the year of data production. The x^-based probability of the asymmetry to be constant is 32% 
and 63% for momentum and time-period dependences, respectively. 

A summary of the systematic errors for the branching ratio ^+ ) is given in Table [3 

The fractional error due to PID and tracking has been estimated as 2.1% of the branching ratio, 
computed as the sum in quadrature of 1.1% for PID and 1.8% for tracking |^. 



Table 4: Summary of systematic errors for the CP asymmetries. 



Source 




(0vr±) 






AcP [10-2] 


AcP [10-2] 


Acp [10-2] 


MC simulation 


0.06 


0.06 


0.06 


Background estimate 


0.63 


0.32 


0.49 


Using r > 6.0 


0.22 


0.15 


0.01 


Using ri > 8.8 


0.46 


0.54 


0.54 


Total 


0.81 


0.65 


0.73 



5 RESULTS and SUMMARY 

In summary, we have searched for a CP asymmetry in — > K~K~^tt~^ ^ — > 4''^'^ i s-iid — > 
K*^K^ decays and measured the branching ratio of K^ir^ decays all with a data sample 
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Figure 2: Mass distributions for positively charged (left) and negatively charge (right) D^g) mesons 
for events satisfying the requirement on r > 4.3. Figures (a) and (b) are for all KKtt candidates, 
while (c) and (d) are for (pn candidates, and (e) and (f) for K K candidates. Figures (g) and 
(h) are for i^vrvr candidates. Signal (yellow or light shaded) and sidebands (red or darker shaded) 
regions are shown for and decays in (a) and (g), respectively. 
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of 79.9 fb collected by the BABAR experiment. 

Preliminary measurements of the CP asymmetries are Acp{K^K^tt^) = (1.4 it 1.0(stat.) it 
l.l(syst.)) X 10-2, Acp{(t>TT^) = (0.2 ± 1.5(stat.) ± 0.8(syst.)) x lO'^, and Acp{K^K*'^) = (0.9 ± 
1.7(stat.) ± 0.8(syst.)) x 10^^. These results are in agreement with previous published results 
with our results in the resonant modes having significantly smaller errors. 

Further, we obtain a preliminary branching ratio for — > K~K~^7r^ decays relative to — > 
decays of (10.7 ± 0.1(stat.) ± 0.2(syst.)) x lO^^. This result is a significant improvement 
over previous measurements 0. 

Table 5: Systematic errors for the branching ratio. 



Source 


Error [10"^] 


PID + tracking 


0.22 


Background estimate 


0.05 


Using r > 6.0 


0.00 


Using ri > 8.8 


0.02 


Total 


0.23 
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